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ABSTRACT 

Quasi-periodic oscillations (QPOs) in the X-ray lightcurves of accreting neutron star and black hole 
binaries have been widely interpreted as being due to standing wave modes in accretion disks. These 
disks are thought to be highly turbulent due to the magnetorotational instability (MRI). We study 
wave excitation by MRI turbulence in the shearing box geometry. We demonstrate that axisymmetric 
sound waves and radial epicyclic motions driven by MRI turbulence give rise to narrow, distinct peaks 
in the temporal power spectrum. Inertial waves, on the other hand, do not give rise to distinct peaks 
which rise significantly above the continuum noise spectrum set by MRI turbulence, even when the 
fluid motions are projected onto the eigenfunctions of the modes. This is a serious problem for QPO 
models based on inertial waves. 

Subject headings: accretion, accretion disks — MHD — turbulence — waves — X-rays: binaries 



1. INTRODUCTION 

Millisecond variability has been observed for some 
time with the Rossi X-Ray Timing Explorer in both 
neutron star and black hole binaries (e.g. McClin- 
tock & Remillard 2004, van der Klis 2004). In ad- 
dition to a broad continuum in the power spectrum, 
QPOs are observed. These QPOs have been modeled 
as hydrodynamic waves in geometrically thin accretion 
disks (for reviews see e.g. Wagoner 1999 and Kato 
2001), or "toru s-like" flows with significant radial pres- 



ZVVi), or torus- like Hows with significant radial pres- 
sure gra£hentej|Giai^ 



Kluzniak et al.ll2004t iLee et al]l2004t iBlaes et.aL 



I200fl) . 

Turbule nce seeded by the magn etorotational instabil- 
ity (MRI; iBalbus fc Hawlevlll99lh is widely believed to 
provide the stresses responsible for angular momentum 
transport and accretion. This raises a series of ques- 
tions. Do waves idealized as hydrodynamic perturbations 
on a laminar background still exist in turbulent, magne- 
tized accretion flows? Does the time-averaged flow act 
as a resonant cavity supporting standing waves? What 
is the steady-state wave amplitude due to turbulent ex- 
citation? The simplest and most controlled geometry to 
explore these quest ions is that of the shearing box (e.g. 
lHawlev et al.l fl995h . One published attempt has been 
made to look for disc rete mode frequenci es in a shearing 
box MRI simulation (Brandenbura 2005), with negative 
results (see figure 7 of that paper). Analytic estimates 
for turbulent exc itation of waves in accretion dis ks have 
been reported in lNowak &; Wagoner! (|1993i fl~995l) . 

We revisit this problem, finding that MRI turbulence 
does indeed excite distinct peaks in the power spectrum, 
and that these peaks can be identified with certain classes 
of hydrodynamic modes, specifically axisymmetric acous- 
tic waves and a global axisymmetric epicyclic oscillation. 
Hydrodynamic inertial wave modes, however, are not 
clearly detected, a serious problem for QPO models that 

1 Kavli Institute for Theoretical Physics, Kohn Hall, University 
of California, Santa Barbara, CA 93106; arras@kitp.ucsb.edu 

2 Department of Physics, Broida Hall, University of California, 
Santa Barbara, CA 93106; blaes@physics.ucsb.edu 

3 Jet Propulsion Laboratory, MS 169-506, California Institute of 
Technology, Pasadena, CA 91109; neal.turner@nasa.jpl.gov 



are based on these modes ( "g-modes" in diskoseismology 
parlance). 

The plan of the paper is as follows. In section [21 we 
review the shearing box and compute the linear wave 
frequencies. We present power spectra from numerical 
simulations in section and discuss the results. We 
present our conclusions and discuss their relevance to 
models for QPOs in section 0] 

2. THE SHEARING BOX 

We solve the magnetohydrodynamic (MHD) equations 
as applied to an isothermal, differentially rotating, mag- 
netized accretion disk, in the unstratificd shearing box 
approximation: 
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Here p, P, v, and B are the mass density, gas pres- 
sure, velocity and magnetic field. The sound speed c s 
is constant. The simulation domain represents a small 
region of the disk near the midplane orbiting with fre- 
quency £1. The non-inertial reference frame of the or- 
bit is taken into account by including a Coriolis force, 
as well as "tidal" forces due to the difference of gravi- 
tational and centrifugal forces, represented by the term 
2qQ 2 xe x , where q = — dlnfl/dhiR is the shear param- 
eter. Cartesian coordinates are used, with the radial, 
azimuthal, and vertical coordinates labeled x, y, z, re- 
spectively. The azimuthal and vertical boundaries are 
periodic, and the radial boundaries are shearing-periodic 
(fluid passing through one radial boundary appears on 
the other at an azimuth that varies in time according to 
the difference in orbital speed across the box). 

An initially weak magnetic field will drive turbulence in 
the shearing box due to the MRI, which is a nearly incom- 
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pressible instability. A second branch of nearly incom- 
pressible perturbations which we call "inertial waves" 
also exists, at least on a laminar background flow. In 
the limit of zero magnetic field, the specific angular mo- 
mentum gradient is the restoring force for these inertial 
waves. To understand the two solutions with a simple 
example, linearize about a time-independent background 
with constant density p and pressure P, constant mag- 
netic field B = B y e y + B z e zi and velocity v = —qQxe y . 
For incompressible, axisymmetric waves with space-time 
de pendence exp(ik x x+ik z z —iu}t), the dispersion relation 
is IjBa.lhus fc HawlevllT99l 
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where the positive sign gives inertial waves and the nega- 
tive sign the MRI. Here «Az — B z / (Airp) 1 / 2 is the vertical 
Alfven speed and k = [2(2 — q)] 1 ' 2 ^! is the fluid epicyclic 
frequency. 4 In the zero field limit, incompressible inertial 
waves have the dispersion relation uj 2 = K 2 k 2 /k 2 . 

Although inertial modes modified by magnetic tension 
exist for the constant background chosen above, once the 
medium becomes turbulent, their possible existence as 
normal modes of oscillation becomes highly questionable. 
Low order inertial waves have frequencies of the same 
order as the energy-bearing MRI "eddies", and hence 
are susceptible to nonlinear interactions with the MRI 
turbulence. 

The simulations described here have zero net magnetic 
flux in each direction. Furthermore, the time and space- 
averaged magnetic pressure is small. We attempt to un- 
derstand waves in such a turbulent background by again 
linearizing about a time-independent background with 
constant density and pressure, velocity v = —qttxe y , 
but now including compressibility and neglecting mag- 
netic fields. In this case, axisymmetric inertial-acoustic 
waves satisfy the dispersion relation 
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For zero wavenumber one recovers the fluid epicyclic 



mode 



In the limit 



<C czk , which is a 



good approximation even for the lowest wavenumbers 
in our simulations, inertial waves again have the disper- 
sion relation co 2 = K 2 k 2 /k 2 , while acoustic waves satisfy 
oj 2 = c 2 k 2 . In a shearing box of dimensions (L x , L y , L z ). 
these waves will occur at discrete frequencies because the 
wavenumbers will be quantized according to 



2im x 



and k? 



2nn 7 
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where n x and n z are integers. 

Because of the shearing boundary conditions, non- 
axisymmetric waves in the shearing box have an 
amplitude and frequency wh ich change with time 
ijGoldreich fc Lvnden-Belil ir965). and therefore do not 
exist as normal modes. The radial wavenumber changes 
with time as k x (t) = k' x + qfl(t — t Q )k' y , where k' x and 

4 Not to be confused with the epicyclic frequency of a particle 
in a gravitational field. 



k' y are the wavenumbers with respect to comoving co- 
ordinates [x, y' = y + qQx(t — to)], and to is the time 
of minimum k x {t). When k x (t) ~ k' x , non-axisymmetric 
waves have similar frequencies to axisymmetric waves. 
For k x (t) 3> k' x , an approximate WKB dispersion rela- 
tion for the time-dependent fr equency is (see eq. 72 in 
IGoldreich fc Lvnden-Bel]|li965j) u(t) ~ c s k' y qVLt, showing 
that the frequency changes on the shear timescale of the 
disk. As k x (t) becomes large, eventually the wave will 
damp away. If the wave has lifetime T, then we expect 
it to give rise to a broad peak in the power spectrum of 
width ~ ck y qQT around the corresponding axisymmetric 
mode. 

3. POWER SPECTRA 

Our simulations use box size (L x , L V) L z ) — (1.25, 4, 1) 
with {N x ,N y ,N z ) = (40,64,32) grid points in each di- 
rection, respectively. The shear parameter q = 3/2, cor- 
responding to Keplerian rotation. The net magnetic flux 
in each direction is zero and the initial weak field has 
a spatial dependence given by B = Bq cos(2irx/L x )e y + 
Bq sm(2irx/ L x )e z . The initial magnetic pressure is 400 
times less than the gas pressure. The sound speed, or- 
bital frequency, and average density have the constant 
values c s = Q = 10~ 3 and p = 1. The fiducial simulation 
is run for 1000 orbits. Averaged over time and the box, 
the total magnetic pressure in the turbulence that results 
is 0.65% of the average gas pressure. However, there is 
significant variance, with small regions of local magnetic 
pressure reaching 10% of the gas pressure. 

Figure n shows "position space" power spectra for den- 
sity, radial velocity, and vertical velocity fluctuations. 
Three spectra are plotted for each. Data recorded at 
a specific point (x, y, z) in the box contain fully non- 
axisymmetric motions. Fluid motions with nodes along 
the j/-direction contain most of the power in MRI tur- 
bulence, hence we also record time series for azimuthally 
averaged data at the same (x, z) and box-averaged data 
in order to increase signal to noise for oscillatory modes. 
We interpolated the data to uniform spacing in time, 
padded the array to a multiple of two, multiplied by 
the Bartlett window function to reduce spectral leak- 
age ijPress et alJll992j) . and used a Discrete Fast Fourier 
Transform to produce the power spectrum. The power 
spectra are normalized to the root-mean-square of the 
time series. To produce figure ^ we binned the power 
into 512 logarithmically spaced frequency bins in the fre- 
quency range shown. This has the effect of multiplying 
the power spectrum by a factor of uj. The frequency uj/2ir 
is expressed in units of inverse orbital time. Arrows show 
the position of low order acoustic and inertial waves us- 
ing eq. 10. While all three panels show curves for non- 
axisymmetric and axisymmetric data, the box-averaged 
power is only visible in the plot for v x , appearing as a 
sharp peak at to = k = f2, the epicyclic frequency. This 
epicyclic mode has been reported i n previous shearing 
box simulations ijHawlev et al.l lT995). Numerical damp- 
ing causes the sharp decline in power at high frequencies. 
In runs at higher resolution, this dropoff was pushed to 
higher frequencies, as expected. 

The most striking feature in figure ^ is the axisymmet- 
ric sound waves, which appear as sharply defined peaks. 
The power in these peaks approaches the broad con- 
tinuum due to non-axisymmetric MRI turbulence, and 
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Fig. 1. — "Position space" power spectrum for density (upper), radial velocity (middle), and vertical velocity fluctuations (lower). In 
each plot, the upper curve is for the non-axisymmetric data and the lower curve is for the axisymmetric data. The sharp peak at the orbital 
frequency in the middle plot is the box averaged data. The frequency is in units of cycles per orbital period. Power has been rebinned into 
logarithmically spaced frequency bins, which effectively multiplies the power spectrum by one power of frequency. Axisymmetric hydrody- 
namic normal mode frequencies from eq. J7J are shown as arrows. The epicyclic mode is at the orbital frequency (unity). Inertial modes 
and sound waves should lie below and above the epicyclic frequency, respectively. From right to left below the epicyclic frequency, arrows 
show the expected axisymmetric inertial mode frequencies for (n x ,n z ) = (1, 2), (1, 1), (2, 1), (3, 1), (4, 1) and (5, 1). From left to right above 
the epicyclic frequency, the axisymmetric acoustic waves have (n x , n z ) = (1, 0), (0, 1), (1, 1), (2, 0), (2, 1), (0, 2), (1, 2), (3, 0), (0, 3), (4, 0), and 
(5,0), respectively. 



clearly rises above the continuum for at least two of the 
peaks. Axisymmetric sound waves with n z = have 
much larger density and radial velocity fluctuations than 
those with n x = 0. This is consistent with the larger con- 
tinuum seen in v x relative to v,. A number of previous 
inves tigations (e.g. IStone et alJll996t IGardiner fc Stone! 
12005(1 have noted non-axisymmetric sound waves in snap- 
shots of the density profile. These non-axisymmetric 
sound waves also have n z = and n x ^ 0, consistent 
with our findings for axisymmetric waves. 

There are no significant peaks observed for inertial 
modes in figure ^ However, this is not proof of their 
non-existence. The epicyclic mode was buried in the 
noise for non-axisymmetric and axisymmetric data, and 
only became visible by box-averaging. Moreover, the ob- 
servable luminosity variation decreases strongly as wave- 
length decreases, due to averaging of hot and cold spots 



on the disk. Hence large lengthscale, small power iner- 
tial modes may in principle be observable as compared 
to short lengthscale, large power MRI eddies. We now 
show that this is in fact not the case. 

Figure [3 shows the "momentum space" power spec- 
trum. Time series for azimuthally averaged v x and v z 
were computed by taking spatial transforms in the x and 
Z directions, for n x and n z nodes, respectively. Hence, 
turbulent eddies and inertial modes are on the same foot- 
ing in this plot, as we compare their power at the same 
lengthscale. The temporal power spectrum for the re- 
sulting time series was computed as in figure ^for 128 
frequency bins in the range shown. The normalization of 
the curves is such that the sum over all n x , n z , and fre- 
quency bins gives the root-mean-square of the time series 
(compare to axisymmetric curves in figure [T]). 

The epicyclic mode is again clearly visible in figure |21 
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Fig. 2. — "Momentum space" power spectrum for axisym- 
metric radial velocity (upper), and vertical velocity (lower) fluc- 
tuations. Each curve is labelled by (n x ,n z ) used for the spatial 
Fourier transform. Power has been rebinned into logarithmically 
spaced frequency bins, which effectively multiplies the power spec- 
trum by one power of frequency. Axisymmetric hydrodynamic 
inertial mode frequencies from eq. J7J are shown as arrows with 
{n x ,n z ) = (0,0), (1,2), (1, 1), (2, 1) from right to left. 



The inertial modes are qualitatively different. Two of the 
lowest order modes show broak peaks of width Suj ~ lu at 
the expected frequencies. Without the analytic mode fre- 
quencies as a guide, it would be hard to recognize these 
peaks as significant. Higher order modes (not plotted 
here) show even less evidence of peaks at the expected 
frequencies. The lack of clear peaks for inertial modes 
suggests that they are strongly affected by MRI turbu- 
lent stresses. We note that although the epicyclic mode 
has frequency similar to inertial modes, it is immune to 
magnetic tension forces as it is nodeless in the shearing 
box, and hence is less susceptible to turbulent buffeting. 

Overall, we find that axisymmetrizing or box-averaging 
the data from simulations greatly enhances the signal 
from oscillation modes, a technique easily adapted to 
global simulations. In addition, long integration times 
will increase the signal to noise of the nearly sinusoidal 
waves relative to the stochastic background from turbu- 
lence. It is not immediately obvious, however, why we 
clearly see peaks in our non-axisymmetric data whereas 
iBrandenburd (|2005f) did not. Perhaps this is because his 
simulation was stratified in the vertical direction. 
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4. CONCLUSIONS 

We find that turbulence driven by the MRI naturally 
excites radial epicyclic motion and axisymmetric sound 
waves, but there is scant evidence of inertial waves. This 
may imply that inertial modes are strongly affected by 
MRI turbulent stresses, and cease to exist as high qual- 
ity oscillators. We note that the zero net magnetic flux 
simulations presented here give rise to the weakest MRI 
turbulence possible. Inertial waves could be even more 
strongly affected in simulations with net magnetic flux. 

Inertial waves, called "(/-modes" in the diskoseismol- 
ogy literature, have been proposed as the origin of some 
obser ved QPOs in X-ray binaries fWagonerl 1 1 999t iKatol 
120011) . However, our results imply that all but the longest 
wavelength of these hydrodynamic modes are destroyed 
by the turbulence. Even the longest wavelength modes 
are substantially altered, perhaps to the point of becom- 
ing unobservable. This presents a serious problem for 
models invoking (7-modes to explain observed QPOs. 

We are not clearly able to pick out non-axisymmetric, 
shearing waves in the analysis done so far. There are 
broadened peaks in the non-axisymmetric data around 
the frequencies of axisymmetric waves, but this is not 
clear confirmation. Furthermore, "real" accretion disks 
may in fact admit standing nonaxisymmetric waves with 
constant pattern speeds. There is substantial global sim- 
ulation data which demon strates this fact, at le ast for 
hydrodynamic disks (e.g. iBlaes & Hawlevl lT988). The 
shearing box cannot capture these modes, because the 
shearing box boundary conditions preclude their exis- 
tence. 

In summary, this paper demonstrates that QPOs can 
be found in simulations of MRI turbulence, at least in 
a shearing box geometry. The well-defined boundary 
conditions of this geometry help ensure the existence of 
modes, but it is not clear that such well-defined bound- 
aries will exist in real accretion disks. Global MRI 
simulations could in principle determine self-consistently 
whether wave cavities exist and contain trapped modes. 
Traveling acoustic wav es have certainly been observ ed in 
such simulations (e.g. iDe Villiers fc Hawlevll2003j) . but 
no evidence for modes at discrete frequencies has yet 
been reported. We encourage extensions of our analy- 
sis techniques to these simulations. 
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